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ABSTRACT 

We present deep radio observations of four nearby dwarf spheroidal (dSph) galaxies, designed to 
detect extended synchrotron emission resulting from weakly interacting massive particle (WIMP) dark 
matter annihilations in their halos. Models by Colafrancesco et al. (CPU07) predict the existence 
of angularly large, smoothly distributed radio halos in such systems, that stem from electron and 
positron annihilation products spiraling in a turbulent magnetic field. We map a total of 40.5 deg 2 
around the Draco, Ursa Major II, Coma Bernices, and Willman 1 dSphs with the GBT at 1.4 GHz to 
detect this annihilation signature, greatly reducing discrete-source confusion using the NVSS catalog. 
We achieve a sensitivity of a su b < 7mJy/beam in our discrete source-subtracted maps, implying 
that the NVSS is highly effective at removing background sources from GBT maps. For Draco we 
obtained approximately concurrent VLA observations to quantify the variability of the discrete source 
background, and find it to have a negligible effect on our results. We construct radial surface brightness 
profiles from each of the subtracted maps, and jackknife the data to quantify the significance of the 
features therein. At the ~ 10' resolution of our observations, foregrounds contribute a standard 
deviation of 1.8mJy/beam < a ast < 5.7mJy/beam to our high-latitude maps, with the emission in 
the Draco and Coma dominated by foregrounds. On the other hand, we find no significant emission in 
the Ursa Major II and Willman 1 fields, and explore the implications of our upper limits for particle 
dark matter using the fiducial models of CPU07. For a WIMP mass M x = 100 GcV annihilating into bb 
final states and B = 1 /iG, the individual dSphs place an upper limit on the annihilation cross-section 
log((crv) x , cm 3 s _1 ) < —25 for one set of charged particle propagation parameters adopted by CPU07; 
this is comparable to that inferred at 7-ray energies from the two-year Fermi-LAT data. We discuss 
three avenues for improving the constraints on (&v) x presented here, which may yield an additional 
order of magnitude in sensitivity for some models. We conclude that deep radio observations of dSphs 
are highly complementary to indirect WIMP searches at higher energies. 
Subject headings: galaxies: dwarf — dark matter — radio continuum: galaxies 



1. INTRODUCTION 

The standard cosmology derived from astronomical ob- 
servations such as the cosmic microwave background, the 
large-scale galaxy distribution, and the kinematics of in- 
dividual galaxies and clusters predicts that the universal 
matter density is dominate d by dynamical l y col d, col- 
lisionless dark matter (e.g. iKomatsu et al.1 feoilf ). Al- 
though there are few clues to its nature, its observed 
abundance requires that a dark matter particle have an 
annihilation cross section (av) X: th ~ 3 x 1 0~ 26 cm 3 s -1 
if it was once in thermal equilibrium (e.g. IPorter et al.l 
1 2 1 If ) . The correspondence of this cross-section to parti- 
cles near the weak scale makes weakly interacting mas- 
sive particles (WIMPs) very attractive dark matter can- 
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didates (se e Uungman et al.l 119961 iBergstroml 120001 and 
lFen^[2010l for reviews). Whether WIMPs constitute the 
dark matter, and whether clues to the WIMP identity 
can be obtained are therefore pressing questions. 

Indirect WIMP searches focus on detecting standard 
model particles that result from WIMP annihilations 
or decays. Given the relative paucity of astrophysical 
sources and the relative robustness of annihilation 
signature predictions at high energies, 7-ray searches 
hav e yielded some of the tightest constraints on (<rv) x 
(see lStrigarill2012l for a recent review). A wide variety 
of objects ha ve been targeted, i nc luding the Galactic 
center (e.g. lAbdo et all I2010al iDobler et al.l 120101: 



Abazaiian fc Hardingl 120121 : lAbazaiian fc Kaplinghatl 
20121 iWenigerl 12013 ). the diffuse Galact i c and extra- 
galactic backgrounds (jAbdo et al.1 l2010d : iCirelli et al.1 
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20"l0f IPapucci fc Strumial l20Tol iBaxter fe Dodelsonl 

2012a) and galaxy cl usters 
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tainties in astrophysical backgrounds, 



However, uncer- 
the dark matter 
distributions of the targeted systems and boost factors 
due to dark matter substruc ture weaken constraints de- 
rived from these studies (e.g.lSu et alj|2010t iDobler et al.l 
l201lt Imouai20ll1: lAckermann et al.l 1 20121 . 

By contrast, nearby dwarf spheroidal galaxies (dSphs) 
are attractive targets for indirect dark matter se arches: 
they are strongly dark matter dominated (e.g. iMateol 
|1998[) . and stringent upper limits on their star forma- 
tion rates imply negligible as trophysical 7-ray emission 
(e.g. lAckermann et al.ll2012bl) . In particular, the line-of- 
sight integrals of the dark matter distributions that are 
consistent with their stellar kinematics (the "J-values") 
predict that the annihilation signal from D raco should be 
the s t rongest among the class i cal dSphs dStrigari et al.l 
12001 ICharbonnieretall 12011 I Walker et al.l 120 111) . J- 
valucs for the ultra-faint dSphs Ursa Major II (hereafter 
UMall) , Coma Berenices (hereafter Coma) and Willman 
1 (hereafter Willi), while more uncertain, may be a fac - 
tor of a few larger (jStrigari et al.l l2008t iStrigarillMTl ). 
The basic optical properties of these four dSphs are given 
in Tabled) 

Accordingly, many 7-ray experiments have tar- 
geted dSphs to search for WIMP annihilations (e.g. 
Abdo et al.ll2010bb lAckermann et aLll2011b [Aleksic et al.1 
20111: lAliu et al.l 120121 ). No detections have been re- 



ported. The strongest upper limits on (av) x for WIMP 
masses M x < 500 GeV stem from observa tions with the 
Large Area Telescope (LAT) on Fermi (|Atwood et al.l 
I2009T ) : analysis of the two-year Fermi-LAT data assum- 
ing M x = 100 GeV and annihilation into bb produces up- 
per limits (ov) x < 10 _25 cm 3 s _1 for individual dSphs, 
and {ov) x < 7 x 10~ 2 6 cm 3 s~ 1 when observa tions of 
10 dSphs are combined (lAckermann et a.l.|[201 lfc see also 
iGeringer-Sameth fc Koushiappasll2011|) . 7-ray searches 
are therefore closing in on (av) Xit h expected for a thermal 
relic at these masses, though sensitivity gains of factors 
of a few are still required. 

A variety of WIMP annihilation channels produce 
non-thermal electrons and positrons that could be de- 
tected when they lose energy through synchrotron, 
inverse Co mpton scattering or bremmstrahlung pro- 
cesses (see IProfumo fc Ullioll2010L for a review). Com- 
pared to 7-ray searches, this approach is compli- 
cated by uncertainties in charged particle propaga- 
tion and energy losses. Nonetheless, because the 
expected signals span the electromagnetic spectrum 
and can be relatively sensitive to (<Jv) x , a variety of 
multi-wavelength WIMP search es have also been car- 
ried in the Galactic center (e.g. IBergstrom et al.l 120061: 
Hooper et all 120071: ICrocker et al.l 120101: iLinden et al l 
20l3P 



20101: iLahaet al.l I2012D. the c 
Hooper! 120081 : iFornengo et al.l 



and galaxy cluste rs (e.g 



background (e.g . 
20121: IHooper et al.1 12012 ) 
Perez-Torres et al.l I200S : 



iColafrancesco et al.l 12001 Istorm et al.l I2012T). In addi 



tion, radio (iTasitsiomi et all 120041: IBorriello et al. | l201( ; 



ISiffert et aljboilD and X-rav (jJeltema fc Profumoll2008 ) 



A series of models bv lColafrancesco et al.l (120071 . here- 
after CPU07) suggest that WIMP annihilations in Draco 
will produce a smoothly distributed, degree-scale radio 
synchrotron halo. Fig.[T]shows the predicted annihilation 
signal at v = 1.4 GHz for a representative set of CPU07 
models with M x = 100 GeV annihilating into bb, a turbu- 
lent magnetic field strength in Draco of B = 1 /iG, and a 
Kolmogorov diffusion coefficient describing charged par- 
ticle propagation of either D — 3 x 10 28 cm 2 s _1 ("set 
#1") or L> = 3 x lO^cn^s- 1 ("set #2"). Because 
the size of the diffusion zone corresponds to twice that 
of Draco's stellar distribution, a crude estimate of the 
predicted annihilation halos for other dSphs is obtained 
by scaling the profiles radially according to their half- 
light radii rn (Table [TJ . This approach assumes a fixed 
annihilation flux, which is reasonable given the similar- 
ity between the J-values of the sy stems considered here 
dStrigari et al.lf2008t IStrigarill2012l) . 

Fig. Q] plots the predicted halo intensities for CPU07's 
optimistic choice of (crv) x ~ 3.4 x 10 _23 cm 3 s _1 , a value 
that is now strongly ruled out by 7-ray searches. With 
this (crv) x , the predicted signals would be easily de- 
tectable with existing single-dish radio telescopes. While 
CPU07's (<rv) x is no longer relevant, their dark matter 
models and particle propagation parameters remain ro- 
bust: realistic values of {av) x can therefore be probed 
with the models of Fig. [T] by scaling them linearly in 
intensity when comparing to observations. 

However, few strong co nstraints on the radio flux den - 
sities of dSphs exist (e.g. iFomalont fc Geldzahlerll 19791) . 
and no dedicated searches for extended radio halos in 
dSphs have been performecQ- The dearth of extended 
emission searches in dSphs stems in part from the ob- 
servational challenges associated with reaching sensitivi- 
ties below the confusion limit of single-dish te lescopes, 
impo sed by discrete background sources (e.g. ICondonl 
119741 ). Radio searches for WIMPs in dSphs therefore re- 
quire a combination of single-dish observations to detect 
the predicted extended emission, and higher resolution 
intcrfcromctric observations to subtract discrete sources 
from the single-dish maps and beat their nominal confu- 
sion limit. 

Accordingly, we have obtained deep radio continuum 
maps of dSphs using the Robert C. Byrd Green Bank 
Telescope (GBT) in order to search for extended radio 
halos resulting from WIMP annihilations. In this paper, 
we present v = 1.4 GHz observations of degree-scale re- 
gions centred on Draco, UMall, Coma and Willi. We 
subtract discrete sources from the Stokes I maps us- 
ing the Karl G . Jansky Very Larg e Array (VLA) Sky 
Survey (NVSS; iCondon et all 119981 ) to reduce discrete- 
source confusion ( £12.11 and H3.ip . and confirm that dis- 
crete source variability is unimportant to our analysis 
with near-simultaneous VLA observations of the Draco 
field ( §2.21 and §3.2| . After accounting for the filtering of 
large-scale flux due to our baselining procedure ( §3.3j) . we 
use the empty UMall and Willi fields to place upper lim- 
its on {<jv) x in the context of the models shown in Fig.Q] 



observations of Local Group dwarf galaxies have been 
used to constrain WIMP properties in this context. 



1 N ote that the limit for Draco from ITfomalont & Gcl dzahlerl 
IH979T ) shown in fig. 12 of CPU07 has little relevance in the con- 
text of dark matter searches: these intcrfcromctric observations 
resolve out structures larger than few arcminutcs and are therefore 
insensitive to degree-scale annihilation halos. 
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Fig. 1. — Predicted radial surface brightness profiles for Draco 
by CPU07, scaled to v = 1.4 GHz and units of mjy per 9.12' beam, 
appropriate for observations with the GBT at this frequency. The 
models assume M x = 100 GeV annihilating into bb with (uv) x = 
3.4x 10 — 23 cm 3 s -1 (now strongly ruled out by 7-ray searches), and 
a Kolmogorov diffusion coefficient of cither Do = 3 X 10 28 cm 2 s — 1 
("set #1") or Do = 3 X 10 26 cm 2 s _1 ("set #2"). The upper radial 
coordinate axis is expressed in units of the Draco half-light radius 
r/j (Table [TJ; we apply this model to other dSphs in our sample by 
scaling it radially to their r^. 

(A color version of this figure is available in the online journal.) 

(SHJ). Wc demonstrate that for a class of dark matter 
models like those of CPU07, deep radio observations are 
highly complementary to 7-ray searches for constraining 
the properties of particle dark matter ( §4.1|) . and describe 
an observational program that could improve the limits 
on (o~v) x found here by an order of magnitude for some 
models ( jjHU) . 

2. OBSERVATIONS AND DATA PROCESSING 

In order to search for extended radio halos in the dSphs 
in Table [TJ we mapped 1.5°- to 4°-square regions cen- 
tered on each dSph with the GBT at a frequency of 
v = 1.4 GHz. This mapping frequency is a compromise 
between synchrotron signal strength and mapping speed 
(which decrease/increase with increasing frequency, re- 
spectively) . This frequency also affords the use of pub- 
licly available NVSS survey data to subtract discrete 
background sources . Be low, we discuss the details of our 
GBT observations fi j2.ip as well as our near-simultaneous 
VLA observations of the Draco field to assess the im- 
pact of discrete source variability in our final GBT maps 

2.1. GBT Observations and Data Processing 

Observations were made with the GBT in numer- 
ous observing sessions under the auspices of pro- 
grams AGBT07C085 (for Draco) and AGBT09A085 (for 
UMall, Coma, and Willi). Maps were made in the on- 
the-fly mode, scanning at 104'/minute, with 2 seconds 
of integration time per pixel, giving a total time for one 
4° x 4° map of 2.8 hours. The scan rate was chosen 
so that with the spectrometer backend integration time 
of 1 second, the telescope motion is 20% of a beam full 
width at half maximum (FWHM) per integration, result- 
ing in negligible beam smearing. RA and dec scans were 
interleaved to reduce the effect of 1// noise, either at- 
mospheric or instrumental, on the final map. The GBT 



Fig. 2. — GBT beam model (solid line) at u = 1.4 GHz used to 
convolve NVSS maps of each field to the GBT resolution in order 
to subtract discrete sources from each dSph field. Above a response 
of / = — 7dB, the GBT beam is well-approximated by a Gaussian 
with a FWHM of 9.12' (dashed line). 

(A color version of this figure is available in the online journal.) 

backends were configured to simultaneously record the 
receiver flux for all polarization products using both the 
Digital Continuum Receiver (DCR) and the GBT Spec- 
trometer. The spectrometer produces both continuum 
and H I information. We defer a discussion of polariza- 
tion and spectral data products to a future publication 
(J. Aguirre et al. 2013, in preparation). In this paper, 
we restrict our analysis to the two DCR XX and YY 
polarization channels, which, when averaged over paral- 
lactic angle (PA), approximate the continuum Stokes I 
intensity. The basic characteristics of the observations 
are presented in Table [5J 

In order to search for extended halos for realistic val- 
ues of (av) x , we need to reach sensitivities below the 
confusion limit of the GBT at v = 1.4 GHz by sub- 
tracting discrete background sources (see fjTJ). The 45"- 
resolution NVSS provides a strong handle on this back- 
ground. Note that although FIRST (Faint Images of the 
Radio Sky at Twenty Centimeters; iBecker et al.l 119951) 
boasts a higher angular resolution and sensitivity per 
beam than the NVSS, it has the significant disadvan- 
tage of resolving out smoothly distributed emission on 
angular scales larger than ~ 1'. We therefore subtract 
only sources detected by the NVSS. The GBT is suffi- 
ciently sensitive that NVSS sources are detected at high 
signal-to-noise in each individual ~ 10-second raster scan 
across a given map. In order to reduce cross-calibration 
uncertainties and to obtain a high-quality subtraction of 
the discrete sources using the NVSS, we directly calibrate 
the GBT data off of the latter. 

For each scan, the raw GBT time-ordered data d(£») = 
di are modeled as 

d i =p i + sNVSS i , (1) 

where Pi is a first-order polynomial (mean and slope) to 
remove baseline drifts in the data, NVSSi is the NVSS 
template (in units of mjy per GBT beam) convolved to 
the GBT resolution, and s is a scale factor determined 
for each scan. Beam-maps constructed from scans across 
bright continuum sources show that the v = 1.4 GHz 
GBT beam is well-approximated by a Gaussian with a 
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TABLE 1 

Basic Properties of the dSph Sample 



Target 
(1) 




(ao,5 ) 
(J2000) 
(2) 


(h,bo) 
(3) 


(arcmin) 
(4) 


Reference 
(5) 


Draco 


17 h 20 m 14.4 s 


±0.6 S , 57°54'54"± 08" 


86.4°, 34.7° 


io.o±°; 3 2 


1 


UMall 


08 h 51 m 29.9 s 


±4.0 S , 63°07'59"± 07" 


152.5°, 37.4° 


14.1 ±0.3 


2 


Coma 


12 h 26 m 59.0 s 


±0.9 S , 23°54'27"± 08" 


241.9°, 83.6° 


5.8 ±0.3 


2 


Willi 


10 h 49 m 21.9 B 


±0.8 S , 51°03'10"± 11" 


158.6°, 56.8° 


9 -3+0.2 


1 



Note. — Col. 2: RA and dec (ooi^o) °f the stellar dSph ccntroid; Col. 3: Galactic co- 
ordinates (ioi^o) computed from (c»o><5o); Col. 4: h alf-light radiu s of an expo nential model 
of the stellar distribution; Col. 5: Reference: (1) IMartin et all ({20081 ): (2) IMunoz et all 

(Unio). 



FWHM of 9.12' above -7dB, but that it falls off more 
rapidly than this at fainter levels. Our beam-maps are 
in good agreement with the GBT ray tracing model of 
Sri Srikant (priv. comm.); we therefore use an azimuthal 
average of this model, shown in Fig. [2j to convolve the 
NVSS data to the GBT resolution. Use of this more ac- 
curate model instead of a Gaussian significantly improves 
point source subtraction. 

The values of pi and s are determined by a least squares 
fit to the given scan's data. We use a first-order polyno- 
mial for pi as a compromise between removing 1 /f noise 
and not removing too much large-scale emission, and 
place quantitative limits on our sensitivity to extended 
halos in ij3.3l The baseline-subtracted, NVSS-calibrated 
data Di are then given by 



D, 



(di -Pi)/s 



(2) 



For conciseness, we refer to maps created from the Di 
as "unsubtracted" because they contain discrete sources. 
By contrast, the baseline-subtracted, NVSS-calibrated 
and discrete source-subtracted data D su b^ (which we will 
call "subtracted") that we require for our analysis are 
given by 



Dsub.i = (d z — pi - sNVSSi)/s 



(3) 



The subtracted scans were visually inspected for dra- 
matic deviations from zero which could result from in- 
terference, atmospheric fluctuations due to poor weather 
conditions, or other problems; a small fraction of scans 
were rejected for these reasons. We also discard scans 
obtained at elevations below 15°, which is the case only 
in the Draco field. The remaining data had XX and YY 
polarizations combined to produce Stokes I, and gridded 
onto the sky using a standard cloud-in-ccll method on a 
40" grid. The maps intended for visual inspection, and 
presented in Figs. [3]-[5j were constructed using a median 
gridding kernel in which the median of the integrations 
within a 3 x 3 pixel box centered on each pixel is adopted 
as the final value. However, all computations were per- 
formed on maps constructed without median gridding, 
to ensure that the map weights were linearly propagated 
throughout the analysis. 

The unsubtracted Stokes I GBT maps obtained from 
Di for the four mapped fields arc shown in the top row of 
Fig. [3J All fields are plotted on the same angular scale, 
and are dominated by discrete sources. This is illustrated 
in Fig. @] for Draco, which shows the inner 2° x 2° of the 
unsubtracted Stokes I GBT map and the NVSS map of 
the same region, convolved to the GBT resolution. The 



striking similarity between the discrete source patterns 
in the two panels suggests that source variability does 
not strongly affect our subtracted maps; we confirm this 
using near-simultaneous VLA observations in i j3.2l 

The second row in Fig. [3] shows the weight maps for 
the Stokes I data, which quantify the relative data con- 
tribution to each pixel on a linear scale from to 1 (black 
and white in Fig. [5j2-[3h, respectively). All calculations 
performed include only map pixels whose weights exceed 
half the median value of the non-zero map weights. We 
verified that our results do not depend on the precise 
value of the median weight threshold adopted. 

The third row in Fig. [3] shows "difference maps" for 
each field. The difference maps are a measure of map- 
ping artifacts, which dominate the uncertainties in our 
analysis (see below). Specifically, we jacknife the data to 
produce three pairs of maps for each field created from 
only half of the collected data: that below ("Split A") 
and above ("Split B") the median time of observation, 
the median telescope elevation and the median telescope 
PA. The morphology of any real sky feature should be 
independent of these observing properties. Accordingly, 
for a perfectly mapped region Split A and Split B for 
each property would be identical, and subtracting Split 
B from Split A would produce a map containing no sky 
signal. We define our difference map as (Split B - Split 
A)/2, so that the variance in these maps has the same 
scaling as that of Stokes I (because halves of the data are 
subtracted to produce the difference map). Note that 
finer data splits (e.g. into thirds) were also attempted, 
but the resulting maps were too noisy to be useful. 

We therefore use the statistics of the features in the 
observed difference maps as a measure of the mapping 
errors. We create difference maps in time, elevation and 
PA for each field. In Fig. [3H3J, we show the difference 
map for each field where the standard deviation of the 
pixels is the largest, and use this map for our error anal- 
ysis in £13.11 For each field, the telescope property jack- 
nifed to produce the plotted difference map is given in the 
bottom-left corner of the corresponding panel in Fig. |31 
PA for Draco, UMall and Coma, and elevation for Willi. 
As expected, the difference maps have more gaps than 
the Stokes I maps of Fig. [3h.-[5tl, because a pixel must 
have data in each jacknifed half to be included in the 
difference map. 

The subtracted Stokes I GBT maps obtained from 
D su b,i for the four mapped fields are shown in Fig. [5J 
and form the basis of our search for extended radio emis- 
sion from the dSphs at each field center. Note that the 
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upper colorscale limit in this figure an order of magni- 
tude smaller than in Fig. [3ji-[3jl: the vast majority of 
the detected flux in the unsubtracted maps steins from 
the discrete source population. Positive or negative fea- 
tures in the subtracted maps represent enhancements or 
depressions relative to the (positive) absolute sky bright- 
ness, which was subtracted during processing by virtue 
of Pi in equation (fTJ) . 

There are no detectable residuals at the locations of 
the bright NVSS sources that were subtracted from the 
maps in Fig. [5l nor is there a correlation between the 
residual map emission and NVSS source locations. As a 
check, we repeat the analysis described in <J3] by blanking 
regions within 1.5 GBT beams of bright NVSS sources: 
our results remain unchanged. Indeed, there are no dis- 
cernible discrete source artifacts in our subtracted maps 
at all. This demonstrates that the NVSS is highly effec- 
tive at reducing source confusion at the GBT resolution. 

We quantify the sensitivity of the unsubtracted and 
subtracted GBT maps of each field in Table [3l As ex- 
pected, the standard deviation a usu b (col. 2) of the pix- 
els in the unsubtracted Stokes I maps well exceeds the 
standard deviation a su b < 7mJy/beam (col. 3) in the 
subtracted maps. In turn, the values of a su b exceed the 
thermal noise expected from the radiometer equation by 
about an order of magnitude, as a result of both mapping 
errors and real sky variations. To disentangle the two, we 
use the statistics of the difference maps in Fig. [3M5J. The 
estimated contribution of mapping uncertainties to the 
subtracted maps is a map (col. 4), and corresponds to the 
standard deviation of the pixels in the difference maps 
of Fig. [3H5J. Assuming that the mapping and sky vari- 
ances are additive, a ast (col. 5) is then the contribution 
to a su b from astrophysical sources. We discuss the inter- 
pretation of a as t as a measure of Galactic foregrounds in 

S21 

The final column in Tablets a measure of the dynamic 
range of the GBT observations, which we define as the 
ratio DR of the peak flux in the unsubtracted maps to 
the standard deviation a su b of the subtracted maps: DR 
is thus an indication of our ability to probe faint struc- 
tures in the field after discrete source subtraction. We 
find DR ~ 140 for the UMall field, and note that the 
lower values for the other fields arise simply because the 
brightest continuum sources therein are fainter. Because 
we find no discrete source artifacts in any of the fields, 
we conclude that DR is only a lower limit to the dy- 
namic range achievable for deep continuum observations 
with the GBT. Our discrete source subtraction is not 
an important contributor to the noise in the subtracted 
maps. It is therefore feasible to produce degree-scale, 
discrete-source subtracted GBT maps with sensitivities 
Csub J; 7mJy/bcam at 1.4 GHz, well below the nominal 
confusion limit. 

2.2. VLA Observations and Data Processing 

To constrain the variability of the discrete sources in 
the Draco field, we obtained VLA observations under 
program AA315 on 2007 November 4, near the middle 
of the 3-month period over which the GBT data for this 
field were taken. The VLA was in B-configuration dur- 
ing that observing cycle, and we therefore configured the 
observing setup to match that of the FIRST survey, with 
one 7-channel, 3 MHz, dual-polarization frequency band 



centered at each of 1365 MHz and 1465 MHz. A compar- 
ison of the measured fluxes from our observations to the 
corresponding FIRST catalog entries provides a conser- 
vative upper limit on the variability in the NVSS sources 
in the same field, which we used for discrete source sub- 
traction of the GBT data (see N3.2p . A total of 5.4 hours 
were spent observing 110 snapshot pointings with the 
VLA in the Draco field, with flux calibrators observed at 
the start, middle and end of the observations and phase 
calibrators observed every 15-30 mins. 

The mapping scheme was designed to reproduce the 
FIRST pointing grid in a 4° x 4° region centered on 
Draco, but an error in the interim VLA mapping soft- 
ware produced smaller pointing offsets than requested in 
RA. The final map therefore spans 17 h 05 m 30 s < a < 
17 h 28 m 30 s and 55°50' < S < 59° 50': this 3° x 4° region 
is centered on Draco in dec, but offset from it by ~ 0.5° 
to the West. The basic characteristics of the observations 
are presented in Table [5J 

The VLA data were reduced in AIPS (|Greisenll2003l) . 
After interactively flagging bad baselines due to inter- 
ference or poor instrumental performance, the flux cal- 
ibrator was used to correct for the bandpass response 
functions of the antennas. The data were then flux and 
phase calibrated using standard AIPS routines, and the 
visibilities at 1365 MHz and 1465 MHz averaged. Each 
pointing was separately imaged and cleaned down to a 
residual noise level of 0.6 mJy, and then corrected for the 
geometric distortion produced by the VLA's snapshot 
mode. The pointings were then mosaicked into a single, 
primary beam-corrected image. 

Source fluxes in the VLA map were measured by fit- 
ting Gaussian components using the AIPS task JMFIT. 
To compare these measured fluxes to the correspond- 
ing FIRST ca t alog e ntries, we follow the procedure of 
Ide Vries et all ([200l : a) we cross-correlate the posi- 
tions of our measured sources with those of the 08Jull6 
FIRST catalog using a matching radius of 3" (roughly 
half a synthesized beam), b) we compare the measured 
peak fluxes of unresolved, matched sources to that of 
their FIRST counterparts, restricting the comparison to 
sourc es brighter than 2mJy (c.f. fig. 2 of Ide Vries et al.l 
l200l . A total of 116 sources meet these criteria: we ex- 
amined each one individually to ensure that none were 
unresolved components of a larger "parent" system. 

Fig. [5] shows the distribution of (measured - FIRST) 
peak flux differences for the 116 sources. The flux differ- 
ences are expressed in units of the net statistical uncer- 
tainty: od = \/ tr„ + Op , where a m is the uncertainty 
returned by JMFIT and ctj? is the RMS map noise at the 
source location listed in the FIRST catalog. Although 
we applied t he FIRST clean bia s correction to our mea- 
sured fluxes (jWhite et al.l 119971 ) , differences in absolute 
calibration and cleaning between the catalogs produced 
a non-zero median flux difference of 0.56 mJy. This has 
been subtracted from Fig. [6] The dotted line in the fig- 
ure shows the best-fitting Gaussian to the distribution. 
The vertical dashed lines denote ±4<t.d, which we adopt 
as our variability threshold. 

We find that 7 unresolved, matched sources in our 
survey area exceed the Aau threshold, and thus exhibit 
variability on the ~ 10-year baseline probed by compar- 
ing our observations to the FIRST catalog; their proper- 
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Fig. 3. — GBT data products. Each column corresponds to a distinct field that is labeled in the top row. All panels span 4° X 4°: the 
horizontal line in the bottom-right corner of a) shows an angular scale of 1°, and applies to all panels. The borders around the Coma 
and Willi fields in the third and fourth columns are blanks around these smaller maps (see Table [2} . Top row: Baseline-subtracted, 
NVSS-calibrated "unsubtracted" Stokes I maps of a) Draco, b) UMall, c) Coma, d) Willi. The linear intensity scale ranges from -10 
to 250 mjy/bcam. Middle row: Weight maps for the Stokes I maps in the row above. The linear colorscale ranges from (black) to 1 
(white). Bottom row: Difference map, created by jacknifing the data, with the largest standard deviation for each field: i), j), k) PA, and 
1) elevation. The linear intensity scale ranges from -10 to 25 mjy/beam. 
(A color version of this figure is available in the online journal.) 



TABLE 2 
Observation and Map Properties 



Field 


Observing Dates 


Intcg.Time 


Map Centre 


Dimensions 


Resolution 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Draco (GBT) 


2007 October - December 3 - 


14.8h 


17 h 20 m , 57°55' 


4° X 4° 


9.12' x 9.12' 


UMall (GBT) 


2009 February - March b 


18. 8h 


8 h 52 m , 63° 08' 


4° X 4° 


9.12' x 9.12' 


Coma (GBT) 


2009 February - March 


8.6 h 


12 h 27 m , 23°54' 


2°. 5 x 2°. 5 


9.12' x 9.12' 


Willi (GBT) 


2009 February c 


1.8 h 


10 h 49 m , 51°03' 


1°.5 x 1°.5 


9.12' x 9.12' 


Draco (VLA) 


2007 November 4 


5.4h 


17 h 18 m , 57°53' 


3° x 4° 


6.8" x 5.3" 



Note. — Col. 1: Field name. Instrument used to obtain the observations described in remaining columns 
is given in parentheses. Col. 2: Dates when majority of data were acquired. Col. 3: Total integration time 
for field. Col. 4: Map centre. Col. 5: Final map dimensions. Col. 6: Angular resolution of final map. 
3 81% of data acquired during 2007 October - December; 19% acquired during make-up sessions in 2009 
and 2010. 

b 91% of data acquired during 2009 February - March; 9% acquired during make-up sessions in 2010 June. 
c 74% of data acquired during 2009 February; 26% acquired during make-up sessions in 2010 June. 
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Fig. 4. — Comparison of the inner 2° X 2° of the Draco field 
from a) our GBT observations and b) the NVSS, convolved to the 
GBT resolution. The horizontal line in the bottom-left corner of a) 
shows an angular scale of 0.5°. In both panels, the linear intensity 
scale ranges from -10 to 250 mjy/beam and the cross denotes the 
stellar centroid of Draco (Table [lj . 

(A color version of this figure is available in the online journal.) 



TABLE 3 
Noise Properties of the GBT maps 



Field 


Cusub 


°~ sub 


C map 


Vast 


DR 




(mjy/bm) 


(mjy/bm) 


(mjy/bm) 


(mjy/bm) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Draco 


33 


6.6 


3.4 


5.7 


88 


UMall 


50 


6.3 


5.1 


3.7 


142 


Coma 


31 


3.6 


1.3 


3.3 


8V 


Willi 


11 


2.3 


1.5 


1.8 


37 



Note. — Col. 1: Field name. Col. 2: standard deviation of 
pixels in unsubtracted map. Col. 3: standard deviation of pixels in 
subtracted map. Col. 4: estimated contribution to cr 3u f, in col. 3 
from mapping uncertainties, or the standard deviation of the pixels 
in the difference map in Fig. [3J-E1- Col. 5: estimated contribution 
to cr su b from astrophysical sources: u^ st = (7^ ub — cr^^. Col. 6: 
Dynamic range: ratio of peak brightness in unsubtracted map and 

CT sub ■ 

ties are given in Table S) The detected variable source 
density is in reason able agreement with the results of 
Ide Vries et all (l200l . who use a similar catalog and ap- 
proach to find 1 variable source per square degree over 
120.2 deg 2 of high-latitude sky on a 7-year baseline: 4/7 
of our variable sources have a fra ctional variabil i ty be - 
low 50% (col. 5 of Tabic 13), while Ide Vries et all (|200l 
report (73 ± 4)%. 
We therefore conclude that on timescales of years, the 




c) Coma 




d) Willi 



Fig. 5. — Discrete-source subtracted Stokes I maps of a) Draco, 
b) UMall, c) Coma, d) Willi. The linear intensity scale ranges 
from -10 to 25 mjy/beam: note that this upper limit is a factor of 
10 smaller than that in Fig,l3hj3tl. The horizontal line in the lower 
left corner of each panel is 0.5° in length, and the cross denotes 
the optical centroid of each dSph (Table [T}. 

(A color version of this figure is available in the online journal.) 
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TABLE 4 

Variable Sources in the Draco Field 



a, <5 


F F 


\F m — 


F F \ 


FR 


(J2000) 


(mJy) 


(mJy) 


(?n) 




(1) 


(2) 


(3) 


(4) 


(5) 


17 h 05 m 04 s .5, 60°03'56" 


3.76 


1.81 


4.53 


1.48 


17 h 19 m 37 s .3, 58°47'55" 


27.3 


7.18 


25.2 


1.36 


17 h 27 m 50 s .7, 57°51'13" 


2.03 


1.40 


5.43 


1.69 


17 h 08 m 02 s .3, 57°44'03" 


4.32 


1.94 


7.22 


1.81 


17 h 06 m 50M, 56°56'15" 


4.33 


2.08 


5.81 


1.92 


17 h 25 m 45 s .7, 56°28'18" 


7.24 


1.22 


4.15 


1.20 


17 h 06 m 57 a .4, 55°54'42" 


3.54 


1.20 


4.08 


1.34 



Note. — Col. 1: FIRST source position. Col. 2: 
FIRST source peak flux. Cols. 3, 4: absolute (measured 
- FIRST) peak flux difference in mJy and in units of the 
net statistical uncertainty. The FIRST clean bias and me- 
dian flux difference for the sample were subtracted. Col. 
5: fractional variability, or the ratio of the brightest to the 
faintest flux (FR > 1). 
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Fig. 6. — Difference between measured peak fluxes of unresolved 
sources in the Draco VLA field and FIRST catalog peak fluxes, 
expressed in units of the net measurement uncertainty crrj. The 
dotted line shows the best-fitting Gaussian to the distribution, and 
the vertical dashed lines denote our variability threshold of Aajj . 
(A color version of this figure is available in the online journal.) 

variable source density and degree of discrete source vari- 
ability in the Draco field is typical of that measured in 
other regions of the high-latitude sky. We investigate the 
impact of this variability on our search for an extended 
radio halo in Draco in 



3. SEARCHING FOR EXTENDED EMISSION 

We wish to assess the likelihood that the subtracted 
Stokes I GBT maps in Fig. [5] contain emission that could 
constitute an extended dark matter annihilation halo. 
To quantify this, we compute the azimuthally averaged 
radial surface brightness profile in each GBT map rela- 
tive to the stellar centroid of the dSph therein, as well 
as the statistical significance of features in the profiles 
fi j3.1[) . We then establish that uncertainties due to vari- 
ability in the discrete source populations subtracted from 
the maps have a negligible impact on the derived profiles 
fi j3.2l) . Finally, we quantify our ability to recover ex- 
tended radio halo profiles given the baselining procedure 
adopted to calibrate the data f £|3.3[) . 

3.1. Radial Profile Derivation 



We carry out the bulk of our analysis on radial profiles 
computed from the subtracted Stokes I maps in Fig. [5j 
The reason for this is twofold. First, because the observa- 
tions are dominated by mapping uncertainties (c.f. i|2.ip 
it is difficult to characterize the noise in each pixel, which 
precludes computing reliable statistics directly from the 
maps. However, uncertainties on radial profile values can 
be robustly determined by jacknifing the data. Second, 
the CPU07 radio halo prediction is azimuthally symmet- 
ric (Fig. [T]), and as a result the comparison of this pre- 
diction to profiles computed from the maps is straight- 
forward. 

We derive radial profiles I(r) from the subtracted 
Stokes I maps by extracting the mean value of all points 
that fall in circular rings as a function of distance from 
the stellar centroid of each dSph (Table [1]). We note that 
the uncertainties on these centroids are much smaller 
than the resolution of our GBT maps, and varying the 
point about which the profile is computed within these 
uncertainties has no impact on the results. 

We exploit the difference maps of each region to ob- 
tain realistic estimates of the uncertainties on each pro- 
file point Ii. We generate 5000 "shifted" realizations of 
each difference map: for each realization, we shift the 
map (and corresponding weights) by a random number of 
pixels in RA and dec, with pixels shifted outside the map 
region wrapped to the other side of it. We then compute 
the radial profile for each realization: for perfect data, 
this radial profile should have I(r) = regardless of the 
structure present in Stokes I. For each difference map, the 
uncertainty on each radial profile value Ii due to map- 
ping errors is then the standard deviation of the points 
obtained from the 5000 realizations. We then adopt the 
largest of the uncertainties on the time, elevation and PA 
difference map profile points as the uncertainty on each 
Ii from the Stokes I maps. 

Fig.[7Jshows radial profiles derived from the subtracted 
maps in Fig. [S] (solid line), as well as those derived 
from the time (dotted line), elevation (dashed line), and 
PA (dash-dotted line) difference maps for each region. 
Qualitatively, all three difference map profiles for each 
field seem consistent with I(r) = 0, as expected if the 
telescope response has been properly calibrated out of 
the data and our error analysis is reliable. The pro- 
files derived from the subtracted maps in the UMall and 
Willi fields (Figs. [7b> and [7b!) also seem to contain lit- 
tle structure. Conversely, the Draco profile (Fig. [TJi) has 
I(r) > (relative to the subtracted mean sky brightness) 
at r < 1.3°, while that for Coma (Fig. [7b) has I(r) < 
at r < 0.7°. 

We quantify the statistics of the radial profiles by 
computing the reduced x 2 statistic for the hypothesis 
that I(r) = 0. We take into account correlated map- 
ping uncertainties by using the covariance matrix derived 
from the profiles of the 5000 difference map realizations. 
Specifically, for each difference map j, the \r statistic is 
given by: 

1 



Xr,Cj 



-1 T T 



N 



(4) 



where Ij is a vector containing N radial profile points 
corresponding to that difference map, Cj is the (N x N) 
covariance matrix computed from the profiles of its 5000 
realizations, and j is one of time, elevation or PA. The 
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value of Xr f° r the Stokes I profile for each field is then 
obtained from: 



^-ICr^iT 

N J 



(5) 



where I is a vector containing N radial profile points 
from one of the Stokes I maps in Fig.O In words, we use 
the difference map covariancc matrix that produces the 
smallest Xr to assess the significance of the features in 
each Stokes I profile. In practice, the Cj used correspond 
to the difference maps shown in Fig. [Sj-O, whose pixel 
standard deviations are largest. 

Table [5] shows Xr obtained for the discrete-source sub- 
tracted Stokes I maps and the difference maps for each 
region. For each entry, the number in parentheses is the 
one-sided p-value of the chi-squarcd test for that entry 
and the number of degrees of freedom in col. 2. Because 
the number of pixels contributing to each profile point 
increases with r, the outermost profile points strongly 
influence the reported statistics. We are most interested 
in quantifying detection statistics near each dSph at the 
centre of the maps, and therefore compute Xr f° r the 
profile points at smaller r than the vertical red arrow in 
each panel of Fig. [7] none of our conclusions change if all 
profile points in each field are included. 

Table [5] confirms that there is no statistically signifi- 
cant structure in any of the difference maps for any field, 
and therefore that our calibration is effective and that 
the difference maps provide a reasonable estimate of the 
mapping uncertainties for each field. We note that the 
statistics of the difference map radial profiles change lit- 
tle if correlated errors are ignored in the analysis: the 
profile points are largely uncorrelated. Among the radial 
profiles for the discrete-source subtracted maps of each 
region, both Draco and Coma show statistically signif- 
icant structure at the 7.7a and 43cr level, respectively. 
As we discuss in Sj4j we attribute this structure to fore- 
grounds. There is no statistically significant structure in 
the UMall field, while that in the Willi field is marginally 
significant at the 5.8cr level. 

3.2. Impact of Discrete-Source Variability 

The similarity between the unsubtracted Stokes I maps 
and the NVSS data for each region (see Fig.|4]for Draco) 
and the absence of residuals at the locations of NVSS 
sources in our subtracted maps (Fig. [5]) suggest that 
discrete-source variability does not strongly impact the 
derived radial profiles for each region. Below, we quan- 
tify this statement using the variability analysis for the 
Draco field from £12.21 

To estimate the importance of variability in the point- 
source population, we conservatively assume that the dis- 
tribution of peak flux differences between our VLA map 
and the FIRST catalog (Fig. [5]) for the Draco field stems 
entirely from intrinsic, low-level variability in the dis- 
crete source population on timescales of years, and that 
this variability is present on the angular scales probed 
by the NVSS. In reality, one expects the variability of 
NVSS sources to be smaller than that of FIRST sources, 
because the extended radio emission resolved out by 
FIRST but retained by the NVSS typically varies on 
longer timescales than probed here. Simulations using 
the variability between our VLA observations and the 



corresponding FIRST map should therefore provide a 
conservative upper limit on the importance of this po- 
tential bias. 

We therefore generate a new realization of the sub- 
tracted Stokes I GBT map for each dSph, and add a 
residual at the location of each NVSS source therein 
whose amplitude is randomly drawn from the gaussian 
in Fig. [6] The probability of drawing a > 4tr outlier from 
this distribution is much lower than required to produce 
the 7 genuinely variable sources in the Draco field (see 
Table @|: to account for this, we also assign the peak 
flux differences of these variable sources to 7 randomly 
selected NVSS source locations in each realization. We 
add these 7 sources to all of our fields, even those that 
are much smaller than the Draco field (Coma and Willi): 
this mimics the effect of a larger space density of vari- 
able sources, which we address below. The ensemble of 
residuals for each realization is convolved to the GBT 
resolution and added to the subtracted map for the cor- 
responding field, and the radial profile is measured. We 
repeat this exercise 5000 times for each field to examine 
the variation in the measured profiles. 

Fig. [5] shows the standard deviation avi at each pro- 
file point computed from the 5000 realizations in each 
of the mapped fields. For all r in all fields, these val- 
ues are at least an order of magnitude smaller than the 
uncertainties on the radial profile points (c.f. Fig. [7]): 
o~vi « o~i. Even in the Coma and Willi fields, where 
the space density of variable sources is assumed to be 
~ 2.6 and ~ 7 times larger than in Draco by virtue of 
the addition of all 7 sources in Table HI the resulting ayi 
increases by at most 30%. It is clear that variability in 
the discrete source population contributes negligibly to 
the uncertainty in the measured radial profiles in Fig. [7J 

The space density of bright (F > 100 mJy) transient 
or strongly vari able radio sources is e stimated to be S, 
0.004 d ee 2 (e.g. IGregorini et all [l98l ICroft et al.l[20Tq 
but see iMatsumura et al.l 120091 ) . and they are therefore 
quite rare in degree-sized fields like those considered here. 
Nonetheless, it is possible for one such event to escape 
detection in our variability study (which only probed one 
~ 10-year baseline) and produce a large residual in the 
discrete-source subtracted map. This scenario seems un- 
likely for our fields, since the largest residuals therein 
does not coincide with the location of NVSS sources. 
Nonetheless, we examined the impact of adding an un- 
resolved residual, with an amplitude equal to the largest 
in the Draco field, at random locations in the subtracted 
Stokes I maps. Unless one is unfortunate enough to have 
this residual fall within ~ 2 GBT beams of the profile 
center, this single large residual does not strongly alter 
the measured radial profiles. In practice, it is straight- 
forward to mask out large residuals when the profile is 
computed, and we therefore consider the effect of this 
uncertainty to be negligible. 

3.3. Sensitivity to Extended Emission 

A key step in processing single-dish radio continuum 
observations is the removal of baseline drifts in the data, 
which we parametrize with pi in equation ([T]). Because 
these drifts have a characteristic scale on the order of 
the map size, an important side-effect of this step is 
that at least some large-scale emission is filtered out. 
This is a well-studied effect of the subtraction of poly- 
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Fig. 7. — Radial profiles measured from the subtracted Stokes I GBT maps in Fig. \E\ (solid line) for a) Draco, b) UMall, c) Coma and d) 
Willi. Note that the horizontal and vertical scales in the panels differ. In each panel, the dotted, dashed and dot-dashed lines show profiles 
derived from the time, elevation and PA difference maps, respectively. For perfectly mapped fields, these profiles would have I(r) = 0. 
Poi nts to the left of the vertical arrow were included in the Xr computation of Tableland the comparisons to the halo models of CPU07 
in £|4,ll see text for details. 

(A color version of this figure is available in the online journal.) 



TABLE 5 
Radial Profile Xr Statistics 



Field 


DOF 


Stokes I 


Time Diff 


El Diff 


PA Diff 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Draco 


9 


7.7(<0.0001) 


1.6(0.11) 


1.1(0.36) 


2.2(0.019) 


UMall 


8 


0.8(0.60) 


2.1(0.032) 


0.5(0.86) 


1.0(0.43) 


Coma 


5 


43(<0.0001) 


1.1(0.36) 


0.1(0.99) 


1.1(0.36) 


Willi 


3 


5.8(0.00078) 


1.6(0.12) 


1.3(0.27) 


1.2(0.31) 



Note. — Col. 1: Field name. Col. 2: Number of degrees of 
freedom in Xr computation. Cols. 3—6: Value of Xr obtained for the 
Stokes I (col. 3), time difference (col. 4), elevation difference (col. 5) 
or PA difference (col. 6) profile. The number in parentheses is the 
one-sided p— value of the chi-squared test given Xr an d the degrees 
of freedom in col. 2. 



nomial or other time-domain functi ons when making 
maps from time ser ies scans (see also lDicker et a,l.ll2009t 
lAguirre et al.ll20"TTl) . We find that fitting a second-order 
polynomial to the data time-series is too aggressive, and 
filters out most of the annihilation halo flux predicted by 
CPU07. We have therefore implemented a linear baselin- 
ing procedure, and evaluate its filtering effect here. 

It is clear that the sensitivity of the final maps to a 
given emission feature is a function of both its character- 



istic scale and morphology: to use the extremes as an ex- 
ample, our maps have full sensitivity to discrete sources 
but are blind to constant emission across the field, as the 
latter would be completely filtered out when the data 
are baselined. The most reliable way to quantify our 
sensitivity to extended halos is therefore to inject them 
directly into the time-ordered data and to attempt to 
recover them post-processing. We carry out this proce- 
dure as follows: we inject each of the profiles in Fig. [TJ 
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Fig. 8. — Estimated uncertainty in the discrete-source subtracted 
Stokes I profiles in Fig. [7] due to variability in the discrete source 
population. The lines show the standard deviation of profile points 
derived from 5000 realizations of the Draco (solid line), UMall 
(dotted line), Coma (dashed line) and Willi (dot-dashed line) fields 
obtained by adding residuals drawn from the Gaussian distribution 
in Fig. \E\ as well as the 7 outliers from Table [4] randomly at the 
locations of NVSS sources. 

(A color version of this figure is available in the online journal.) 

scaled radially to the r/, of the dSph in each field (Ta- 
blc[T]), directly into the time-ordered data di to produce a 
"data+halo" input di + hi. We then process di + hi in the 
same manner as for di in i j2.ll The maps produced from 
di are subtracted from those produced from di + hf. to 
zeroth order, this eliminates thermal and mapping noise 
as well as the input di to produce an image of hi after 
processing. Finally, we compute radial profiles from the 
processed hi maps in the same manner as for the di for 
each field. 

We verified that the reconstructed surface brightness 
of the processed halo images is linear with input halo 
surface brightness to within 2% for 5mJy/beam < I < 
100 mJy/beam, as expected given the success of our dis- 
crete source subtraction [c.f. equation. (JTJ]. We therefore 
produced processed images for a single "set #1" and "set 
#2" halo predicted by CPU07 for the dSph in each field. 

Figs. I§1 and ITU1 illustrates the impact of our processing 
pipeline on different input halos. The solid lines therein 
show input "set #1" and "set #2" halo profiles, and the 
dotted lines show the corresponding radial profile com- 
puted from the processed halo maps. As expected, our 
sensitivity to a given profile depends on its radial extent 
relative to the map size: our processing pipeline recov- 
ers 100% of the "set #2" halo flux predicted for Draco 
(Fig. [10]) . but only ~ 60% of the peak surface brightness 
for "set #1" . Because a larger fraction of the "set #2" 
flux is recovered, we are correspondingly more sensitive 
to (<jv) x computed from this profile than for "set #1": 
our maps will therefore yield a range of constraints on 
(o~v) x by virtue of the propagation parameters adopted. 

Note that in Figs. I§1 and ITU1 the processed "set #1" 
halo profiles are negative at large r: this is another ar- 
tifact of the baseline subtraction. By removing a mean 
and slope from the maps, the linear baselining routine ef- 
fectively imparts a "curvature bias" relative to the map 
center. Note, however, that this curvature bias does not 
imply a positive bias in large scale flux. The curvature 
bias could explain why the non-zero radial profiles for 
Draco and Coma, which we attribute to foregrounds be- 



low, appear to peak roughly at the map center. We ex- 
plore mapping strategics to mitigate this and the filtering 
effect described above in £14.4.21 

Throughout this paper, we compare all data products 
to the processed halo profiles (c.f. the dotted lines in 
Figs. IH1 and [TU|) to account for any filtering or curva- 
ture bias due to our baselining procedure when placing 
constraints on the extended emission from the targeted 
dSphs. 

4. DISCUSSION: IMPLICATIONS FOR WIMPS AND 
FUTURE PROSPECTS 

With measured radial profiles f q3.1[) , the impact of 
variability explored f ^3.2p and an assessment of our sen- 
sitivity to extended radio emission in hand f ^3.3[) . we 
turn now to the implications of the subtracted Stokes I 
maps for WIMPs. Specifically, we interpret our derived 
radial profiles in the context of the CPU07 predictions in 
fc|4.11 we address the limitations imposed by foregrounds 
in i )4.21 and discuss the plausibility of invoking large- 
scale magnetic fields in dSphs in fc|4.3l In i j4.4[ we exploit 
the lessons learned in this study to suggest an observing 
program that should improve our sensitivity to {o~v) x by 
an order of magnitude for some WIMP models. 

4.1. Comparison to CPU 07 Models 

We restrict ourselves to comparisons between our data 
and the CPU07 predictions shown in Fig. [TJ scaled ra- 
dially to each dSph according to (see JT). While 
these assumptions are a gross over-simplification of the 
detailed modelling carried out by CPU07, they nonethe- 
less provide a first order estimate of the (o~v) x that are 
probed by our maps. 

The radial profiles in Fig. [7] and statistics of Table [5] 
show that there is no significant structure in the UMall 
field, while that in Willi is only marginally significant. 
These dSphs thus lie in regions with low astrophysical 
foregrounds even at our sensitivity and resolution, and 
exhibit no evidence for extended radio emission. Assum- 
ing that the emission from the dSphs in these fields is 
consistent with zero, we use the errorbars on the profile 
points in each field to place upper limits on the strengths 
of the radio halos allowed by the data. 

Fig. IH1 shows the result of this exercise. The hatched 
regions in each panel show the uncertainties on our ra- 
dial profile measurements, and the dotted lines in each 
panel show the CPU07 "set #1" and "set #2" predic- 
tions - processed as described in H3.3I - normalized to 
be consistent with the profile uncertainties at 95% confi- 
dence. Because the outermost profile points pull strongly 
on the fit but are of little astrophysical interest, we in- 
clude only points to the left of the vertical arrow each 
panel in the fit, as in the Xr computations of Table [5j 
The solid lines show the corresponding input halo pre- 
dictions, from which we infer the upper limits on {o~v) x 
given on the right-hand side of the panel. 

Fig. |9] illustrates that our upper limits on (o~v) x from 
UMall and Willi differ by almost two orders of magni- 
tude between the CPU07 "set #1" and "set #2". This 
arises because the latter is more centrally concentrated 
than the former, which has two effects. First, for a given 
total flux, a centrally concentrated profile is easier to de- 
tect than one distributed over more pixels. Second, our 
baselining algorithm filters out less flux for a centrally 



12 



Spckkcns ct al. 



concentrated profile relative to a more extended one (see 

gap . 

The radial profile for Draco shows a statistically sig- 
nificant feature with the same basic morphology as that 
predicted by CPU07. However, it is premature to asso- 
ciate this emission with dark matter annihilation. First, 
the subtracted Stokes I map in Fig. [SJi exhibits a hori- 
zontal band of emission roughly centred on Draco, which 
is suggestive of foregrounds. Indeed, the amplitude of 
the emission seen in the Draco field is within a factor of 
two of the depression in the Coma field center, suggesting 
that the features in both maps stem from astrophysical 
foregrounds on similar angular scales. Recall that the 
correlation of these features with the map center stems 
at least in part from the curvature bias imparted by our 
baselining procedure ( §3.3j) . 

Finally, there is tension between the {o~v) x implied by 
the Draco profile if foregrounds are negligible and the 
upper limits that we derive from the UMall and Willi 
fields. Fig. [POl shows the best-fitting processed "set #1" 
and "set #2" predictions to the Draco profile points, the 
corresponding input predictions and the implied {o~v) x . 
It is clear that if one associates the majority of the emis- 
sion in the Draco field with an extended halo, the {o~v) x 
required exceeds our upper limits from UMall and Willi. 
We therefore conclude that the Draco field is likely con- 
taminated by foreground emission. Because of the small 
sample size probed here, we defer a thorough, simultane- 
ous analysis of the predicted annihilation halos and fore- 
ground contributions to a future publication (see £ 14.2 .2[) , 
and do not attempt to constrain {o~v) x from either the 
Draco or Coma fields here. 

Fig. [9] shows that for models of the type predicted by 
CPU07, J(r) - UMall and Willi fields imply (av) x < 
10- 23 cm 3 s- 1 for "set #1" and (av) x < lCT^cm^- 1 
for "set #2" at 95% confidence. For "set #2" , our upper 
limit on {<jv) x is commensurate with the two-year Fermi- 
LAT limits for individual dSphs (Ackermann et al1l2011l : 
iGeringer-Sameth fc Koushiappasl 1201 ID . and within ~ 
0.5 dex of the thermal prediction \og{av) t x = —25.5 (see 
ijpj). We caution that a direct comparison between the 
Fermi-LAT (o~v) x and our limits is muddled by differ- 
ences in halo properties and WIMP models adopted, as 
well as uncertainties in the magnetic field and charged 
particle propagation parameters inherent in our con- 
straint. Nonetheless, the comparisons of Fig. [9] suggest 
that deep radio observations of dSphs have the potential 
to provide interesting constraints on particle dark matter 
properties. 

4.2. The Impact of Foregrounds 

In this section, we focus on the limitations imposed by 
foregrounds on radio constraints of (o~v) x , and invert the 
problem in our maps to produce high-resolution, high- 
sensitivity measures of Galactic foreground emission at 
1.4 GHz ( 34.2. ip . We then validate our field selection 
for obtaining limits on {o~v) x in §4.11 using a Bayesian 
estimate of the halo and foreground signals implied by 
the data (j £j~2~2|) . 

4.2.1. Mitigating Foregrounds 

A potential avenue for mitigating Galactic foregrounds 
- which predominantly stem from synchrotron emission 



- is to explore the spectral dependence of the features 
in the subtracted maps of Fig. [5j The spectral index 
of the annihilation signature depends on the adopted 
propagation parameters (c.f. fig. 12 of CPU07), but it 
generally resembles a ~ —0.7 expected for Galactic syn- 
chrotron emission for 1 GHz < v < 10GHz, and steepens 
to a < -1 for 10 GHz < v < 100GHz. So in principle, 
deep v > 10 GHz radio observations would enable spec- 
tral index analyses that may help to identify foregrounds 
and subtract them from the maps. In practice, however, 
mapping degree-scale fields at such high radio frequencies 
to the required sensitivity is prohibitively expensive. 

We did attempt to examine extant survey data 
for the fields in Fig. [5] at lower frequencies from 
the compilation of diffuse G alactic radio emission by 
Ide Oliveira-Costa et all (|2008| ). Their predicted all-sky 
map at v = 1.4 GHz is largely an extrapolation of 
the lHaslam et all (|1982ft v = 408 MHz survey data, 
which has an angular resolution of 1° and an aver- 
age zero-level uncertainty of 3 K. Processing these maps 
through our pipeline as described in £)3.3I we find fluc- 
tuations on_rtie_order of lOm Jy/beam. However, scal- 
ing the lHaslam et al.1 (|1982l ) zero-point uncertainty to 
;/ = 1.4 GHz assuming synchrotron foregrounds, these 
fluctuations lie well within the noise. There are therefore 
no extant foreground maps that have either the sensitiv- 
ity or resolution to compare to our subtracted maps. 

Another possibility for mitigating foregrounds is to 
use the Galactic H I spectral line signals that we 
collected simultaneously with the continuum measure- 
ments presented here (see M2.lt J. Aguirre et al. 2013, 
in preparation). H I is highly correlated with dust, 
which motivates the use of GBT-obtained H I maps 
to clean maps of Galactic foreg round dust emission 
(| Planc k Collaborat ion et al.1 120111 ) There is a correla- 
tion between the synchrotron and dust emission for en- 
tire galaxies which persists on very large scales, but the 
correlation on small scales is expe cted to be tight onl y 
for the highest radio frequencies dBennett et al.l 12003ft . 
Thus we do not expect that the H I emission will be a 
complete representation of foregrounds in our maps. 

We can therefore turn the problem around and as- 
sume that there is no extended radio emission from any 
of the targeted dSphs; then, the maps in Fig. [5] pro- 
vide a constraint on high-latitude Galactic foregrounds 
at unprecedented sensitivity and resolution. In par- 
ticular, a ast in Table [3] is a measure of Galactic fore- 
grounds at 10' resolution and 35° < b < 84°. We find 
1.8 mjy/beam < cr ast < 5.7 mjy/bcam for our fields with 
the largest a ast at the lowest b as expected. However, 
the correlation between a as t and b is not perfect, likely 
due to the small number of fields mapped. Nonetheless, 
our measurements represent an improvement of a factor 
of ~ 50 in sensitivity and ~ 40 in angular resolution 
over current foregro und measurements at 1.4 GHz, such 
as that computed by Ide Oliveira-Costa et al.l (|2008f) . 

4.2.2. A Bayesian Halo and Foreground Estimate 

The Draco and Coma fields were discarded from the 
analysis in §4.11 under the hypothesis that I(r) ^ 
therein stemmed from foreground contamination rather 
than extended radio emission. Although the foreground 
argument seems to be the most plausible interpretation 
of the available data, it raises the uncomfortable possi- 
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bility that I(r) ~ in the UMall and Willi fields results 
from a fortuitous cancellation of a positive halo signal 
and negative foreground fluctuations: in this case our 
limits on {crv) x have little meaning. We regard this sce- 
nario as unlikely but attempt here to simultaneously an- 
alyze all four mapped fields - allowing for a contribution 
to I(r) from both a CPU07 annihilation halo and a fore- 
ground contribution - in order to further validate the 
approach in ^4.11 

We adopt a Bayesian approach using the best-fitting, 
processed CPU07 predictions to each Stokes I profile (as 
in Fig.[TU]for Draco) and their uncertainties. We assume 
that the best-fitting intensities are drawn from a normal 
distribution with the following components: i) a positive 
definite annihilation halo signal, with a common normal- 
ization for all dSphs, ii) an RMS foreground signal which 
can be either positive or negative, and iii) our estimated 
zero-mean measurement noise. For the foreground signal 
we consider two cases: a) a constant foreground across all 
fields, and b) a foreground with a 1/ sin(6) dependence. 
We then find the most likely values of the halo normal- 
ization and the RMS foreground contribution given the 
best-fitting processed "set #1" and "set #2" profiles to 
the measured Stokes I profiles in Fig. [7] for all four fields 
jointly. 

We find that the analysis strongly rules out the zero 
foreground case: it is not possible to interpret the Stokes 
I profiles from all four fields in the context of the CPU07 
models without invoking a foreground contribution. This 
is not unexpected given 7(r) < in the Coma field 
(Fig- Eh); but the analysis nonetheless lends some quanti- 
tative support to the approach adopted in H4.ll We find 
that an RMS foreground amplitude of 4—60 mJy/beam is 
allowed by the data at 95% confidence, with a most likely 
value of 6— 9 mJy/beam. This estimate is consistent with 
1.8mJy/beam < a ast < 5.7mJy/beam obtained in §2. 11 
though with considerably larger uncertainty because an 
extended halo is simultaneously modeled and because the 
Bayesian estimate assumes only one foreground fluctua- 
tion in each map (on the scale of the annihilation halo 
signal), whereas a as t is computed over several indepen- 
dent beams. The most likely value of the halo normal- 
ization peaks near zero for both "set #1" and "set #2", 
justifying our computation of upper limits on (crv) x from 
the UMall and Willi profiles. However, including a fore- 
ground contribution weakens the upper limits on (crv) x 
by ~ 0.5 dex relative to those derived in §4.11 for both 
foreground models considered. 

While the results of our Bayesian analysis are encour- 
aging, we find that our sample of four dSphs is too small 
to meaningfully distinguish foregrounds from the signal 
of interest in a simultaneous fit: our approach in §4.11 
of considering only fields that are unlikely to be con- 
taminated by foregrounds is therefore preferable for the 
present study. Moreover, we have not considered the po- 
tential spatial dependence of the foreground signal within 
the maps, which is poorly constrained at their sensitivity 
and resolution ( §4.2.1[) . Combined with the shortcomings 
of adopting a single halo normalization for all four fields, 
we find that a comprehensive Bayesian analysis is beyond 
the scope of this paper, but is likely to be both feasible 
and preferable as both larger samples and better models 
become available (see ^4.4|) . 
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Fig. 9. — Upper limits on {av) x from a) UMall and b) Willi 
for the CPU07 models shown in Fig. [T] scaled radially using 
for that dSph (Table [l}. In both panels, the hatched region shows 
the radial profile uncertainties from Fig. [JJ The dotted black (red 
in the online version) and grey (green in the online version) lines 
show the CPU07 halo profile shapes for "set #1" and "set #2" 
respectively, processed through our pipeline to account for filtering 
by our baselining procedure and scaled in intensity to be consistent 
with the profile uncertainties at 95% confidence. Only uncertainties 
on profile points to the left of the vertical arrow are included in 
the fit. The solid lines in the same color show these "set #1" 
and "set #2" predictions before processing through our pipeline, 
and \og(av) x given on the right-hand side of the panels is the 
corresponding limit implied by the CPU07 models for a thermal 
relic, in cm 3 s _1 . 

(A color version of this figure is available in the online version of 
the Journal.) 



4.3. Magnetic Fields in dSphs 

Unlike their counterparts at 7-ray energies, WIMP an- 
nihilations in dSphs are detectable in the radio only if 
these systems harbor large-scale magnetic fields. CPU07 
note that the dependence of the signal on the magnetic 
field strength is not trivial because it affects both the 
particle propagation and the synchrotron emission, but 
the sensitivity of their fiducial models scales roughly as 
(<jv) x oc B (see their fig. 14). The magnetic field strength 
therefore has an important influence on the detectability 
of the predicted signal. Here we discuss the plausibility 
that the dSphs considered here have turbulent magnetic 
fields with strengths B ~ 1 [iG, as assumed in the fiducial 
CPU07 models of Fig. ffl 

Direct m easurements of B are difficult to obtain (see 
iBeckl 120111 . for a review) . A standard technique is to 
search for polarized radio emission attributable to syn- 
chrotron radiation, and to use the equipartition theorem 
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Fig. 10. — Comparison between the radial profile computed for 
the Draco field and predictions from the CPU07 models shown 
in Fig. [T] The solid line with error bars shows the radial surface 
brightness profile points as in Fig. [7] The dotted black (red in the 
online version) and grey (green in the online version) lines show the 
best fitting CPU07 halo profiles, processed through our pipeline to 
account for filtering by our baselining procedure, to the data for 
"set #1" and "set #2", respectively. Only profile points to the left 
of the vertical arrow are included in the fit. The solid lines in the 
same color show these "set #1" and "set #2" predictions before 
processing through our pipeline, and log(av) x given on the right- 
hand side of the panels is the corresponding value implied by the 
CPU07 model, in em's" 1 . Note that the detected emission likely 
stems from foregrounds and not WIMP annihilations; see text for 
details. 

(A color version of this figure is available in the online journal.) 

to estimate the field strength. Under modest assump- 
tions, measurement of the synchrotron spectral index 
and the degree of polarization allow one to infer the to- 
tal magnetic field strength. Using this approach, turbu- 
lent magnetic field strengths of ~ 10 /j,G have been found 
in some actively s tar-forming, present-day dwarf irregu- 
lar (dlrr) g alaxies (jChvzv et al.ll2000L [2003 : iKeplev et all 
[20101 l20ll while those i n mor e quiescent s ystems like 
the L MC (|Gaensler et all 120051) and SM C (IMao et ah 
2008 ) are in the ~ 3 — 5 fiG range (see also iChyzv et al 



201 lh . However, the lack of a detectable ISM precludes 



using this approach in prcscnt-day dSphs, and B therein 
is largely unconstrained. 

In principle, an estimate of the line-of-sight integral of 
the magnetic fields and electron densities in dSphs can 
be obtained by measuring the Faraday rotation of the 
polarization angle of the polarized emission from back- 
ground galaxies along the line-of-sight, i.e., the rotation 
measure. This measurement is challenging since most 
radio galaxies have polarized fractions of ~ 2%, and we 
are looking for a change in rotation measure of order 
a few. This requires well-calibrated measurements at 
widely separated frequencies to determine the phase rota- 
tion. Very low frequency telescopes can measure multiple 
wrappings within their bandwidth, but lack the requisite 
sensitivity at present (e.g., LOFAR, PAPER, MWA). It 
is therefore unlikely that Faraday rotation studies will 
constrain dSph magnetic fields in the near future. 

Notwithstanding the difficulties in measuring B in 
dSphs today, evidence for similar evolutionary histo- 
ries between dlrrs and dSphs imply that the latter had 
magnetic fields in the past. The broadband photomet- 
ric properties and resolved stellar populations of nearby 
dwarfs suggest that dlrrs and gas-poor dSp hs had similar 
star formation histories until a few Gyr ago ([Calura et al.l 



l2008t iWeisz et al.ll201ll) . Numerical simulations indicate 
that "tidal stirring" , a combination of tidal and ram- 
pressure stripping by the Milky Way halo, i s effective at 
transforming dlrrs into (classical) dSphs (see lMaverll2010l 
for a review), again suggesting a common origin for these 
two classes. If dSphs and dlrrs share a common origin, 
then dSphs may once have had turbulent magnetic fields 
of order 5 — 10 /zG as well. How this early magnetic field 
evolved after star formation ceased and their gas was 
stripped remains an open question. 

It is nonetheless plausible that a tenuous, ionised ISM 
that has insofar escaped detection the dSphs sustains 
a weak magnetic field. For example, Draco could con- 
tain interstellar gas with a mass that is ~ 10% of that 
in its stars and still sat isfy available upper b ounds on 
its diffuse Ha intensity ([Gallagher et al.1 120031 ). On the 
other hand, if we invoke cquipartitition in UMall and 
Willi and make standard assumptions regarding the syn- 
chrotron path length (I = r/j from Table [1} and the 
proton-to-electron ratio (K a = 100), the upper limits on 
the peak "set #1" and "set #2" flux de nsities in Fig ~ 
would allow B ~ 1 fj,G in these systems (jBeck fc Kraus 
120051 ) . If the radio flux from low- luminosity galaxies is 
suppressed relative to that in L* galaxies (e.g. Bell 2003), 
B could be larger. We conclude that although detailed 
modeling of the magnetic field evolution in dSphs is re- 
quired to place quantitative limits on their most likely 
present-day strengths, it is plausible that these systems 
harbor a turbulent magnetic field with B = 1 /j,G as in- 
voked by CPU07. 

4.4. The Way Forward 

Comparisons between the upper limits implied by the 
UMall and Willi surface brightness profiles and the "set 
#2" CPU07 models in Fig. Q] require that log(cn;) x < 
—25 for each system, which rivals the constraints ob- 
tained by Fermi-LAT for individual dSphs ( ij4.ip . It is 
therefore clear that deep radio observations of dSphs are 
highly complementary to WIMP searches at higher ener- 
gies. We now discuss potential avenues for improving the 
constraints on (o~v) x from radio searches of dSphs rela- 
tive to the results presented here: carrying out a joint 
analysis of a large sample fi j4.4.ip , mapping larger, offset 
fields ( H4.4.2p , and exploring a variety of WIMP annihi- 
lation models ($4X3}. 



4.4.1. A Large Sample of dSphs 



In analogy 



to 



jo int dSph an alyses at 

7-ray energies (lAckermann et ail 120111 ; 

iGeringer-Sameth fc Koushiappasl 120111 ), we antici- 
pate that combining subtracted Stokes I maps for 
several dSph fields should improve our constraints on 
{av) x . The sensitivity of our subtracted maps is limited 
by both mapping errors and Galactic foregrounds (cr map 
and a as t in Tabled respectively): particularly because 
of the latter issue, it is unlikely that deeper observations 
of a single target will improve upon the (o~v) x limits 
presented here. However, these uncertainties should de- 
crease when a larger sample of dSphs is jointly modeled. 
Our exploratory Bayesian analysis in H4.2.2I implies that 
simultaneous modelling of extended annihilation halos 
and the Galactic foreground is feasible with a large 
sample. Indeed, all but one of the 10 dSphs considered 
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by lAckermann et all ()2011f ) are accessible to the GBT, 
as arc several other Galactic dSph s from the recent 
compilation of iMcConnachid (|2~012). It is therefore 
feasible to increase the dSph sample size by a factor of a 
few relative to that presented here to improve the limits 
on (av) x . 

4.4.2. Mapping Larger, Offset dSph Fields 

While the constraints on (<jv) x from CPU07's "set #2" 
models and our data are competitive with those from the 
two-year Fermi-LAT data, the limits from "set #1" are 
considerably weaker ( §4.1|) . As explained in 33.31 this is 
due in part to the filtering effect of our baselining proce- 
dure on predicted halos with sizes comparable to the map 
size. This effect can be mitigated by mapping larger re- 
gions around each dSph. While costly in observing time 
(c.f. Tabled]), obtaining larger maps for several dSphs 
is nonetheless feasible with the GBT. As also discussed 
in §3.3[ it is advantageous to offset the fields relative to 
the dSph centroids in order to minimize the curvature 
about the map center introduced by our baselining pro- 
cedure. We note that larger offset maps of Draco and 
Coma may circumvent some of the limitations imposed 
by foregrounds in those fields, because the baselining pro- 
cedure would be both less sensitive to foreground features 
in the map and less apt to correlate those features with 
the dSph centroids. 

4.4.3. Detailed dSph and WIMP Models 

As emphasized in 34. 1[ our comparisons to the CPU07 
models in Fig. [1] are sufficient to estimate the magnitude 
of (<jv) x probed by our maps. However, our deep radio 
observations clearly afford a more sophisticated analysis 
than the preliminary work presented here. The applica- 
tion of a detailed suite of models fo r each dSph, similar 
to tho se presented by CPU07 and Ueltema fc Profumol 
(|2008| ). is required to thoroughly probe the region in 
M x — {av) x — B parameter space that is consistent with 
the subtracted Stokes I maps for a variety of WIMP 
annihilation channels. 

Considering the three avenues discussed above and the 
lessons learned from this study, we anticipate an improve- 
ment in our upper limits on {<jv) x by a factor of ~ 0.5 dex 
by carrying out a joint analysis of > 10 dSphs that are 
offset by ~ 0.5° from their respective map centers. If map 
sizes are selected to minimize the filtering of flux from 
halos with extents comparable to CPU07's "set #1" , we 
should gain an additional factor of 2 ~ 0.3 dex in sen- 
sitivity to {cv) x for those models. Combined with the 
application of a more sophisticated suite of predictions 
for each dSph, it might be possible to gain up to an or- 
der of magnitude of sensitivity to (crv) x for models that 
resemble CPU07's "set #1". For at least some plausi- 
ble combinations of M x , B and charged particle prop- 
agation parameters, we would therefore begin to probe 
\og(av) x , t h ~ —25.5 expected for a thermal relic. We are 
actively pursuing all of these observational and theoret- 
ical avenues. 

5. CONCLUSIONS 

Models by CPU07 predict that for plausible values of 
the turbulent magnetic field strength and charged parti- 
cle propagation parameters in the Draco dSph, WIMP 



annihilations in its dark matter halo will produce a 
degree-scale synchrotron radio halo that is accessible to 
current singlc-dish facilities. We present deep GBT ob- 
servations at 1.4 GHz of a total of 40. 5° 2 around the 
Draco, UMall, Coma, and Willi dSphs to detect this 
annihilation signature. 

We search for extended radio emission associated with 
the dSphs at sensitivities <j su b J$ 7mJy/bcam ( §2.1[) by 
subtracting discrete sources in the Stokes I maps using 
the NVSS catalog. The subtracted map noise is not 
limited by our discrete source subtraction; it is there- 
fore possible to map well below the nominal confusion 
limit with the GBT when the NVSS is used to remove 
background sources. We obtained near-concurrent ob- 
servations of the Draco field with the VLA ( 32.2p . and 
use them to demonstrate that variability in the discrete 
source population has a negligible impact on our results 

For each subtracted Stokes I map, we compute radial 
profiles about the dSph at the map center, jacknifing the 
data in observation time, telescope elevation and tele- 
scope PA to determine robust uncertainties on each pro- 
file point ffl3~l~j). We find that I(r) ^ with high sig- 
nificance in the Draco field and Coma fields. While the 
radial profile computed from the Draco field resembles 
that expected for an extended halo, several lines of ev- 
idence suggest that the emission stems at least in part 
from foregrounds ( 34.ip . Indeed, our subtracted maps 
probe foregrounds at unprecedented sensitivity and res- 
olution ( 34. 2ft : we find that the standard deviation in 
our maps attributable to foregrounds is 1.8mJy/beam < 
Cast < 5.7 mJy /beam, where a ast roughly anti-correlates 
with the galactic latitude range 35° < b < 84° of the 
maps. 

We find no statistically significant features in the 
UMall and Willi profiles. We use the uncertainties on 
the profile points in those fields to place limits on the 
annihilation cross-section {crv) x in the context of the 
CPU07 models for Draco ( 34.ip . which we scale radially 
according to the stellar distribution sizes of the other 
dSphs. We compare the fiducial CPU07 models with 
M x = 100 GeV annihilating into bb final states, B = 1 /j,G 
and two sets of charged particle propagation parame- 
ters, and account for the filtering of those models by our 
baselining procedure ( §3.3|) . 

In the context of the CPU07 models, our up- 
per limits from the UMall and Will fields imply 
log((cri;) x , cm 3 s _1 ) < —23 for the "set #1" propagation 
parameters, and \og((av) x , cm 3 s _1 ) < —25 for "set #2" 
( 34. ip . The "set #2" constraint is comparable to the lim- 
its for individual dSphs obtained at 7-ray energies from 
the two-year Fermi-LAT data. We discuss three potential 
avenues for further improving the limits on {crv) x from 
deep radio observations: carrying out a joint analysis of 
a larger sample of dSphs, mapping larger fields that are 
offset from the dSph optical centroids, and applying a 
more sophisticated suite of models to the data. For at 
least some of the models considered here, we may gain 
up to an order of magnitude in sensitivity by pursuing 
these avenues. We therefore conclude that deep radio ob- 
servations are highly complementary to indirect WIMP 
searches at higher energies, and have the potential to 
probe (crv) Xt th expected for a thermal relic. 
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